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Abstract. We report a simulation study for bottle-brush polymers grafted on a rigid 
backbone. Using a standard coarse-grained bead-spring model extensive molecular 
dynamics simulations for such macromolecules under good solvent conditions are 
performed. We consider a broad range of parameters and present numerical results 
for the monomer density profile, density of the untethered ends of the grafted flexible 
backbones and the correlation function describing the range that neighboring grafted 
bottle-brushes are affected by the presence of the others due to the excluded volume 
interactions. The end beads of the flexible backbones of the grafted bottle-brushes do 
not access the region close to the rigid backbone due to the presence of the side chains 
of the grafted bottle-brush polymers, which stretch further the chains in the radial 
directions. Although a number of different correlation lengths exist as a result of the 
complex structure of these macromolecules, their properties can be tuned with high 
accuracy in good solvents. Moreover, qualitative differences with "typical" bottle- 
brushes are discussed. Our results provide a first approach to characterizing such 
complex macromolecules with a standard bead spring model. 
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1. Introduction 

Macromolecules with comb-like architecture, where linear or branched side chains are 
grafted onto a backbone chain have found much interest in recent years due to their 
physical and biochemical properties, which offer important benefits (TJ |2} El HI El El [7] . 
Although their structure is complicated, such macromolecules are being synthesized 
rather successfully [H [21 HJ El E]- Atom transfer radical polymerization has been 
efficiently applied to the synthesis of various bottle-brush brushes (also known as 
comb-on-comb brushes) [21 HI El [El [11], where bottle-brush molecules are grafted 
onto a linear chain (or a point (single monomer) in the case of star bottle-brush 
brushes) which serves as the backbone of the macromolecule. Brush polymers have 
been studied extensively as far as it concerns experiments [U [21 El IU El El H21 [HI [H] 
as well as theoretical modeling [ISl US! HZl HH [IHl Efll EH E21 [23]. Already, the 
study of "typical" bottle-brush polymers (flexible linear polymeric chains grafted onto 
a backbone) has attracted attention for potential applications [H E]. For example, 
some work was motivated by the use of these cylindrical brushes as building blocks in 
functional supramolecular structures; applications for actuators and sensors have been 
also discussed, since the structure of these stimuli-responsive polymers can change when 
external parameters such as pH of the solution, temperature etc. are varied [2^1 [25] . 
The study of such effects has already been the subject of previous simulation studies 
of bottle-brush macromolecules with one or two types of side chains where a rigid 
backbone has been considered as a starting point for the study of these complex 
macromolecules [261 [271 [IU [291 ETil EH [321 [331 EH EH [36], where it has also been 
shown that the radial density profiles for bottle-brush polymers with stiff and flexible 
backbones are similar when good solvent conditions are assumed. Bottle-brush brushes 
could serve as a candidate for applications already suggested for bottle-brush polymers, 
but in the case of bottle-brush brush polymers the higher number of varied structural 
parameters could help tuning the resulting properties of such macromolecules with 
a higher accuracy and in various ways. Also, towards understanding more complex 
structures of macromolecules which exist in biological systems, simulation models would 
allow for a better understanding of their complex structure-properties relation [371 EH] • 
Bottle-brush brushes exhibit rich structural behaviour providing even ideas for new 
applications due to the higher multitude of the correlation lengths which is expected 
from their structure. More interesting formations could be expected under poor solvent 
conditions as for typical bottle-brushes [23 [27J [23 1211 EQ] , but for bottle-brush brushes 
more difficulties, such as equilibrating such a systems, could be expected. Moreover, 
the interlay of the local conformation of side chains and the global configuration of 
the backbone adds more complexity to their intricate behaviour. The interpretation of 
experimental data on bottle-brush brushes is expected to be controversial as is already 
shown in the case of bottle-brushes with flexible linear side chains [321 E2j. Also, 
simulations has shown that the application of scaling concepts in the latter case is 
proven problematic even in the case where only a single type of side chains occurs, and 
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good solvent conditions are assumed [321 E3J EH ESJ EE]- Moreover, the application 
of theoretical arguments for the description of these macromolecules has been also 
proven very difficult [TJ EJ [5j [6]. A comparison of the results of computer simulations 
could rather be realised on the basis of effective exponents, as for "typical" bottle 
brushes [31]. In view of the experimental interest for bottle-brush brushes and related 
complex macromolecules some fundamental understanding on their static properties 
can be achieved rather well by computer simulations. In this work we describe large- 
scale molecular dynamics simulations of an off-lattice model of a bottle-brush brush 
under good solvent conditions to provide a first approach to understanding their overall 
structural properties. 

The outline of this paper is as follows. In Sec. 2 we describe the simulation 
model. Then, Section 3 discusses the analyzed properties and our numerical results, 
while Section 4 presents and summarizes our conclusions. 



2. Computational details and simulation method 

To simulate a bottle-brush brush polymer, we use a molecular dynamics method, where 
the monomers are coupled to a heat bath [39|. All the monomeric units are modeled by 
the standard bead-spring model [40j HU |42j |43] , where all monomers are treated as beads 
of mass m. This model has been extensively used in previous simulations of brush-like 
systems and its detail discussion has been given elsewhere [4T]. Here, we give a brief 
description of the parameters we use in our simulations. For the interaction between 
any two monomers a truncated and shifted Lennard- Jones potential acts given by 

\^ LJ [{^r-C—f}+C for r<r c 
U LJ (r) — \ T T (1) 

I for r > r c 

where r c = 2 1 / 6 <tlj is the cut-off of the potential, and the constant C is defined 
such that Ulj{t = r c ) is continuous at the cut-off. Henceforth units are chosen such 
that eu = 1, o~lj = 1, &b = 1, and the mass m of the beads is also taken as unity. The 
connectivity of the beads is guaranteed by the "finitely extensible nonlinear elastic" 
(FENE) potential: 

^FENE(r) = -hr 2 Q ln[l - (-) 2 ] < r < r , (2) 
2 r 

where r = 1.5, k = 30 and Ufene{j") = oo outside the range written in equation (J2J). 
The equation of motion for each bead reads 

d d^L _ , > , „ , 

m ^ = - w --^ + r -W' < 3 > 

In this equation 7 is the bead friction, T^t) describes the random force of the heat bath 
and Ui is the potential each bead experiences due to the presence of the other beads, 
when they are below the cut-off distance. The random forces satisfy the standard 
fluctuation-dissipation relation 

< TM.Tjit) >= 6fc B T 7 ^5(t - t) (4) 
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Figure 1: (Colour online) Schematic representation of bottle-brush brush polymer and 
parameters describing the structure of such system, a is the grafted density that the 
bottle brush polymers are grafted onto the rigid backbone, 7V& is the length of the flexible 
backbone. Also, flexible side chains of length N s are grafted onto the flexible backbone 
with a grafting density o~ s . 

where T is the temperature and ks is the Boltzmann constant. Following previous work 
[39l[40l[4ll[42l[43l[44l[45l|46] 7 = 0.5 and T = 1.2. Here, r = {m LJ a 2 LJ / ' e L j) xl<1 is the 
natural time unit, with units that have been given above. We use the the molecular 
dynamics package LAMMPS [37] to simulate our systems where the equations of motion 
for each bead (equation [3]) are integrated with the velocity-Verlet scheme [3H] with a 
time step At = 0.008. Periodic boundary conditions in the x-direction are applied, 
which is the axis of the rigid backbone where the brush chains with a flexible backbone 
are grafted regularly with a distance 1/er between them. In this study we have considered 
different grafting densities, i.e., a = 0.25,0.50, and 1.00, which corresponds to grafting 
every bead, every second bead, and every forth bead of the rigid backbone. Smaller 
densities for our range of chain lengths would suppress any effects due to the density, 
while higher grafting densities than the ones we consider here would impose difficulties 
in equilibrating our systems. In y— and z— directions, periodic boundary conditions are 
applied as well, but the considered linear dimensions of the simulation box were chosen 
large enough, so that never any interaction of the grafted bottle brush polymers with 
their periodic images could occur. We have considered a variety of parameters as they 
are schematically described in figure [U The backbone where the chains are grafted 
with a grafting density a is rigid, while the backbones of the grafted brushes onto this 
rigid backbone are flexible as well as their side chains. The number of monomers of 
the rigid backbone was N = 128. We have also performed simulations for other choices 
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Figure 2: (Colour online) Two characteristic snapshots of bottle-brush brush polymers 
with different set of parameters. These snapshots show the different structures of bottle- 
brush brush polymers under good solvent conditions. This is similar to the behaviour 
of typical bottle-brush polymers. In the left snapshot the set of parameters a = 0.25, 
Nb = 12, N s = 3 and o~ s = 0.50 was chosen and for the right snapshot, a = 1.00, 
Nb = 48, N s = 24, a s = 1.00. Our choice of parameters show the different structures 
that may result ranging from low densed brushes to homogeneous cylinders. 

of backbone lengths in order to check the influence of periodic boundary conditions on 
the resulting properties. We found that our results are not affected by the presence of 
periodic boundary conditions for the rigid backbone length that is taken here. However, 
if one tries to simulate such a system under poor solvent conditions, then one would 
expect an influence of the periodic boundary conditions on the properties of the grafted 
brushes [26j |27j |28j EHJ [30] . In this manuscript, we will present only results for iV = 128. 
The backbone of the bottle-brushes attached onto the rigid backbone is iV& = 12, 24, 36, 
and 48 (the part of the chains denoted with red color in figure [1]). Higher lengths are 
prohibitely difficult to study. Also, it would be difficult for experiments to access such 
lengths as it has been discussed already for typical brushes [32]. For the side chains of 
the grafted bottle-brushes onto the rigid backbone we have typically considered chain 
lengths shorter than iVj,, i.e., N s = 3,6, 12, and 24. These shorter side chains (denoted 
with blue color in figure [T|) are grafted onto the flexible backbone of the brushes with 
a grafting density a s = 0.5 or 1.0. As shorter chain lengths are used for N s , we did 
not study the grafting density a s = 0.25. Figure [2] shows (left part) typical snapshots 
for the case of a = 0.25, iV& = 48, N s = 3 and a s = 0.5. In the right part the case of 
er = 1.00, Nj, = 48, N s = 24 and o~ s = 0.5 is shown. For our choice of parameters we 
observe structures ranging from low densed brushes to homogeneous cylinders. 

We note here that equilibration of our systems required running our simulations for 
every brush for a time range of 20 x 10 6 MD time-steps. Then we collect 2000 samples 
running the simulations 4 x 10 6 MD steps further. By studying correlations functions 
for the structural properties of the brushes that vary slowly with time, we find out that 
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Figure 3: (Colour online) Density profile p(r) as a function of r perpendicular to the 
rigid backbone for bottle-brush brush polymers with iV& = 12, N s = 3 and cr s = 0.5, for 
the grafting densities a = 0.25,0.5 and 1.0. r is the distance from the rigid backbone. 
Inset shows results for the same set of parameters, but a s — 1.0 



this effort was enough to obtain reliable results. In the following, we present our results 
discussing overall structural properties for these complex macromolecules providing an 
insight for the overall behavior of these polymers. 



3. Results and Discussion 



First, we focus on the effect of the grafting density a on the density profile in the radial 
directions (y — z plane with the rigid backbone extending along the x direction) for 
moderate values of Nj, and N s . Figure [3] shows the density profiles of the system in 
radial directions (y, z) for the case JV& = 12, N s = 3 and different grafted densities 
a. From figure [3l the increase of the grafting density a shows that the chains overally 
become more stretched in the radial direction as it normally happens for typical bottle- 
brushes. For low grafting density a, the density curves decay almost exponentially for 
a = 0.25 and their height reflects analogously the corresponding grafting densities (the 
second peak from the left). Also, by increasing the grafting density a the latter extension 
in density becomes higher as the side chains of the grafted brushes get stiffened due to 
the excluded volume interactions. Moreover, for a = 1.0, figure |3] shows that the density 
profile persists over longer distances from the rigid backbone and indicates that the side 
chains of the grafted bottle-brush molecules tend to stretch in the radial directions. In 
addition, the increase of number of side chains grafted onto the flexible backbone a s 
(inset) hint the effect that the curves of density profile extend further, but the overall 
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Figure 4: (Colour online) Variation of the density profile p(r) as a function of r for the 
case of cr = 1.00, iV& = 12 and four different values of N s (a s = 0.5). Inset shows the 
corresponding results for a = 0.25. 

behaviour of the systems does not change. Also, for the case a = 1 and a s = 1 a 
third peak is formed due to the layering effect as the density increases close to the rigid 
backbone. 

Stretching of the flexible grafted bottle-brushes in the radial directions is more 
pronounced when one changes the length of the side chains of the grafted bottle-brushes. 
Figure H] shows differences in the density profile which result from the variation of N s . 
The increase of N s strengthens the extension of the density profile in the radial direction. 
The density now stays almost constant at rather high distances r (the density remains 
constant even at longer distances from the rigid backbone) and the resulting formation 
of the bottle-brush brush resembles a homogeneous cylinder as it is that of figure [2] 
(right part). For cr = 1 we see this rather constant density up to a distance 7 from the 
rigid backbone (disregarding of course the first two peaks, which are reminiscent of the 
layering effect observed for fluid particles in a box close to the wall). As N s increases, 
the density profile extends to higher distances from the rigid backbone. This extension 
becomes a rather small effect as the grafting density a is low enough that the beads can 
not fill the space close to the rigid backbone under good solvent conditions where the 
chains stretch in the radial direction (inset). 

The variation of the density profile as a function of the radial distance r with iVj, 
is shown in figure [5] (a). This typical graph is for the case of grafting density o = 0.25 
and N a = 24. This "low" grafting density allows us to discuss the particular effects 
avoiding the structures where the density close to the rigid backbone is high. At small 
distances r from the rigid backbone, the difference in the density profiles can not be 
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Figure 5: (Colour online) Density profile p(r) plotted versus the radial distance r. (a) 
The case a = 0.25, N s = 24 and four different brush chain lengths iV& is shown, (b) 
Same as figure HI for the case a = 0.25, Nt> = 48 and different values of N s . a s = 0.5 in 
both cases (a) and (b). 



seen. However, at higher distances, the effect of long side chains (N s = 24) shows higher 
density at larger distances r. For iV& < 24 the density profile looks as it is rather expected 
even for a typical bottle-brush polymer. But for iV& > 24 (and N s = 24) the density 
profile shows a completely different behaviour when the length of the side chains (N s ) 
is also rather high. This behavior is characteristic for the bottle-brush brush polymers. 
The density far from the rigid backbone is higher compared to a typical bottle-brush 
polymer. Without the side chains of the flexible brush, we would have expected and 
almost exponential decay in the density profile. However, as the length N s increases 
(figure HJ) , the density in the outer region (far from the rigid backbone) of the bottle- 
brush brush polymer increases, and this difference is more pronounced, when the length 
Nf, is high enough in order to be far from the almost constant density region close to 
the rigid backbone (figure [5] (b)). 

More interesting is to study the behaviour of the ungrafted end beads of the bottle- 
brushes which are grafted on the rigid backbone in order to check how close they can 
come to the rigid backbone. Figure |6] shows our results for the distribution of these 
end beads and the dependence with grafting density a. Since a increases the peak of 
Pe{t) becomes sharper and thinner. This indicates that the end beads of the flexible 
backbone are moving within a narrower space in the radial direction and they rather 
never reach the region close to the rigid backbone. This effect clearly shows that the 
chains are stiffened due to the increase of the grafting density a and as a consequence, 
the peak is higher in the case of high a. The stretching of the brush chains in the radial 
direction is also seen here and corroborates the results of figure [3j According to figure El 
as a increases, the peaks are positioned at higher distances r. This effect is similar 
to the case of typical bottle-brush polymers where the increase of the grafting density 
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Figure 6: (Colour online) Density profiles of the untethered end beads of the flexible 
backbones (bottle brush chains) p_g(r) as a function of radial distance r from the rigid 
backbone for different values of a. This typical illustration is for a system with jVj, = 48, 
N s = 24 and o~ s = 0.5 (a s = 1.0 in inset). 



induces stretching of the side chains in the radial directions and resulting in the increase 
of the zone which is not accessible to the end beads. The increase of a s strengthens the 
aforementioned effects, i.e., the brush chains become more stiffened and more stretched, 
although the chains are under very good solvent conditions. 

The corresponding results for the end beads distribution for figure i.e., figure [7J 
shows that the increase of Nf, shifts the peaks of the resulting density profiles to the 
right, while these peaks become less sharp, lower and smoother, showing that the flexible 
backbone becomes less mobile in the radial directions even for this case where the length 
of the side chains of the brush N s is rather high. Looking at figures [6] and [7J we can see 
that all the changes of parameters have a rather proportional effect on the properties. 
This points out that carefull tuning of the structure of the bottle-brush brushes can 
result from the design of our macromolecules, which can be performed in a variety of 
ways, due to the variety of the parameters that describe their structure. 

When o~ s increases, the differences in the curves are higher and one can see that the 
curves of the density profiles corresponding to different Nf, do not overlap any more. The 
latter effect points us out that the increase of the grafting density a s has a significant 
effect and it is related analogously to the probability for the flexible backbone untethered 
end of coming close to the rigid backbone. Figure El confirms this conclusion. When 
o~ s = 0.5, for different values of N s , we can almost observe the same height of the peaks 
in the density profiles for the end beads, the curves become slightly narrower, but the 
differences are not as pronounced as in the case of a s = 1.00. This clearly shows that 
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Figure 7: (Colour online )Same as figure [61 variation of density profile of the untethered 
end beads versus radial distance r for the case a = 0.50, N s = 24, and a s = 0.5 (cr s = 1.0 
as inset), and different values of iV&. 
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Figure 8: (Colour online) (a) Same as figure El the flexible backbone (brush chains) 
ends Pe(j) vs. r for a = 1.00, N b = 48, and a s = 0.5 ((b) a s = 1.0), and different values 
of N s = 3(D), 6(0), 12(A) and 24(y). 
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Figure 9: (Colour online) The orientation correlation function G(5x) as a function of 
the correlation distance Sx along hte rigid backbone for a bottle-brush brush polymer 
with parameters 7V& = 12, N s = 3 and er s = 0.5 (cr s = 1.0 in inset) for different grafting 
density o. Only a range of Sx is shown. 



the change in the grafting density a s has a significant effect on the bottle-brush brush 
polymers, which we could not distinguish in the overall radial density profiles. We 
should also note that the increase of a s introduces noticeable differences in the height 
of the peaks, showing that this parameter can play an important role in the resulting 
behaviour of our system showing that the variations in N s result in more pronounced 
differences in the peaks for high grafting density of side chains (a s = 1.0), and they now 
have a better correspondence to the values of N s . 

Figures [9] and [JO] show the correlation functions of orientations of different bottle- 
brushes in order to monitor the extent over which correlations between the different 
grafted brushes on the rigid backbone persist [27]. Such a property has been mainly 
considered for bottlebrush polymers with two types of monomers [27] in order to 
characterize the extent over which a Janus structure exists. However, in our case it 
would provide us the with information of the extent that these compact objects (bottle- 
brushes grafted on the rigid backbone at high densities) are affected by the presence of 
their neighboring bottle-brushes. For each chain the vector pointing from the grafting 
site (the monomer on the rigid backbone where the brush is grafted) to the center of 
mass of the respective grafted bottle-brush is considered. Projecting this vector into the 
2/z-plane and defining a unit vector Si for the i-th bottle-brush, we define the correlation 
function 

G(5x) = (S t -S,). (5) 
As we can see from figure El the influence of the grafting density a for values 0.25-1.00 
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Figure 10: (Colour online) (a) The correlation function G(5x) as a function of the 
distance 5x along the rigid backbone for a bottle-brush brush polymer with fixed grafted 
density a = 0.25, Nb = 24 and a s = 0.5 for different side chain length N s (inset, 
a s = 1.0). (b) G(5x) as a function of Sx for a bottle-brush brush polymer with fixed 
grafted density a = 0.25, N s = 3 and a s = 0.5 with different chain lengths iVj, (inset, 
a s = 1.0). The first peak is at \jo. 

is small, and only some difference can be observed when the value of a changes from 
0.25 to 0.50. Further increase of the grafting density a seems to have a small effect on 
the chains. These conclusions are also confirmed for other set of parameter iV& and N s 
(not shown here). We should emphasize that these first peaks occur at a distance 1/cr as 
expected. Furthermore, the same behaviour is observed when the grafting density a s is 
increased (inset in figure [9]). What we see is that the neighbour bottle-brushes influence 
mainly their first neighbours. Two neighbouring brushes grafted on the rigid backbone 
tend to orient in different directions due to the excluded volume interactions, whereas 
the second neighbours, the third, etc., are almost unaffected from the presence of the 
other grafted bottle-brushes on the rigid backbone. Only in the case that a = 1.00 some 
correlation is slightly detected also for the second neighbors. Deviations for G(Sx) from 
zero for third, forth, fifth, etc. neighboring grafted brushes (figure [9]) are within the 
statistical error. Of course, this is due to the short chain lengths we have considered 
for our flexible brushes. When the radius of gyration of the grafted brushes would be 
higher, these effects play more important role. 

For the set of parameters of figure [101 we have varied parameters Nj, and iV s . Here we 
show only results for the case a = 0.25, in order to discuss the particular effects. Similar 
conclusions are drawn for the other cases. Also the increase of N s (figure ITDTa)) shows 
that the correlations increase, but for N s > 6 the occurring differences are smaller. The 
variation of Nb (figure [TUl (b)) leads to the same conclusion. For the results of figures [9] 
and [TOj we can further observe that the increase of o~ s results in the increase of the 
correlations (which of course takes a negative value for neighboring brushes) along the 
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rigid backbone, as the density of the monomers increases. The increase of a s from 0.50 
to 1.00 increases significantly the density of the monomers between the backbone of the 
grafted bottle-brushes, effect which is seen in the results of this correlation function. But 
even in the case of long lengths N s and iV b and high grafting densities grafted bottle- 
brushes further apart form the first neighbors are hardly influenced. To conclude, for the 
the small lengths 7V& and N s considered here, the main effect comes from the variation 
of the grafting density. For other combinations of parameters (not shown here in order 
to save space) the conclusions remain the same. 

4. Summary 

In the present study, we performed molecular dynamic simulations of a standard bead- 
spring model of flexible bottle-brush polymers grafted onto a rigid backbone under good 
solvent conditions. Our investigation was based on the analysis of overall properties of 
interest, i.e., the density profiles in the radial directions (which is the perpendicular 
direction to the rigid backbone where the bottle-brush polymers are grafted with one 
of the ends of their flexible backbone). We discussed that the increase of any of the 
parameters a, a s , 7V&, and N s results in the stretching of the bottle-brushes in the radial 
directions although we are under good solvent conditions, as has been also discussed 
for typical bottle-brushes with a rigid or a flexible backbone. We can also clearly see 
by analysis of the correlation function G(Sx) that mainly only first neighbor grafted 
bottle-brush chains are affected due to the excluded volume interactions between their 
monomers, even for our extreme cases where the grafted densities a and o~ s , and the 
lengths Nj, and N s are high. As the size of the grafted brushes increases, this effect 
is expected to be more pronounced. The strength of the excluded volume effects 
between two neighboring bottle-brush polymers grafted onto the rigid backbone can be 
tuned in various ways, due to the high number of structural parameters characterizing 
such macromolecules, compared to a typical bottle-brush polymer under good solvent 
conditions, where the only varying parameter (considering also a rigid backbone) is the 
length of the grafted linear polymeric chains and the grafting density a. We found 
that the variation of these parameters affect in a proportional way the properties of 
the macromolecules. That is a change of one parameter for all the grafted brushes has 
an analogous measurable effect to the structural properties of these macromolecules, 
although at a first glance their structure seems complicated to be studied by computer 
simulations. Thus, the properties of these macromolecules can be tuned efficiently by 
the change of parameters, of course when these are in a good solvent. The lengths of 
the chains in our study were rather short, but we should stress here that such lengths 
could be also accessible in the synthesis of bottle-brush brushes. We also find that for 
the range of parameters studied here, the untethered ends of the chains have not any 
possibility of being anywhere close to the rigid backbone even for small a. The distance 
and the range that these end beads (possible carriers of some particular substance) can 
access close to the backbone can also efficiently be tuned, which is rather impossible 
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to control with such accuracy in a typical bottle-brush with flexible side chains. The 
width of this zone (unaccessible to the end beads) depends strongly on the N s and a s 
and on the increase of the density depending on the other parameters. We hope that our 
results will stimulate further study of these macromolecules and the work of analytical 
predictions. 
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